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Destruction of a solid rocket stage of a launch vehicle can create a thermal radiation 
hazard for an aborting crew module. This hazard was assessed for the Constellation 
Program (Cx) crew and launch vehicle concept. For this concept, if an abort was initiated in 
first stage flight, the Crew Module (CM) will separate and be pulled away from the 
malfunctioning launch vehicle via a Launch Abort System (LAS). Having aborted the 
mission, the launch vehicle will likely be destroyed via a Flight Termination System (FTS) in 
order to prevent it from errantly traversing back over land and posing a risk to the public. 
The resulting launch vehicle debris field, composed primarily of first stage solid propellant, 
poses a threat to the CM. The harsh radiative thermal environment, caused by surrounding 
burning propellant debris, may lead to CM parachute failure. A methodology, detailed 
herein, has been developed to address this concern and to quantify the risk of first stage 
propellant debris leading to the thermal demise of the CM parachutes. Utilizing basic 
thermal radiation principles, a software program was developed to calculate parachute 
temperature as a function of time for a given abort trajectory and debris piece trajectory set. 
Two test cases, considered worst case aborts with regard to launch vehicle debris 
environments, were analyzed using the simulation: an abort declared at Mach 1 and an 
abort declared at maximum dynamic pressure (Max Q). For both cases, the resulting 
temperature profiles indicated that thermal limits for the parachutes were not exceeded. 
However, short duration close encounters by single debris pieces did have a significant effect 
on parachute temperature. Therefore while these two test cases did not indicate exceedance 
of thermal limits, in order to quantify the risk of parachute failure due to radiative effects 
from the abort environment, a more thorough probability-based analysis using the 
methodology demonstrated herein must be performed. 


Nomenclature 

A t = Area of Nylon 

A 2 = Area of spherical debris 

C = Specific heat of nylon 

CM = Crew Module 

Cx = Constellation 

D = Magnitude of D 

D = Vector of distance between parachutes and a debris piece 

F 12 = View factor 

FS = First Stage 

FTS = Flight Termination System 

LAS = Launch Abort System 

LOC = Loss Of Crew 

LOX = Liquid Oxygen 

Max Q = Maximum dynamic pressure of a trajectory 

MET = Mission Elapsed Time 

Q = Energy 

q = Heat flux 

r = Radius of a spherical debris piece 
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T = Temperature 

e = Emissivity of parachutes 

0 = Stefan-Boltzmann Constant 

p N = Mass to Area ratio of the parachute nylon. 

01 = Angle between normal vector of A t and the vector D 

0 2 = Angle between normal vector of A 2 and the vector D 

I. Introduction 

A solid rocket stage on a launch vehicle can create a thermal radiation hazard for an aborting crew module. This 
hazard was assessed for Constellation’s (Cx) Ares 1 launch vehicle and Orion crew vehicle. The propulsion system 
of Ares I consists of a first stage solid rocket motor and a liquid oxygen (LOX) and liquid hydrogen upper stage. If 
an abort is initiated during first stage flight, the Orion Launch Abort System (LAS) is fired to quickly pull the Orion 
Crew Module (CM) away from the launch vehicle in an effort to save the crew. During this abort there is an option 
to detonate the "headless” launch vehicle for public safety. Detonation of the Ares I results in thousands of 
fragments thrown outward from the destruct location. In such a scenario, there is a risk that launch vehicle debris 
could strike the CM and lead to a Loss Of Crew (LOC) event. The majority of these debris fragments originate from 
the first stage propulsion system and are comprised of burning solid rocket propellant. The risk of debris strikes in 
first stage abort scenarios have been studied in depth by the Cx Program. However, a secondary debris concern 
exists - the thermal effects of burning solid rocket propellant on the CM, and specifically on the parachute system. 

An abort trajectory analysis on Ares I was performed using a first stage debris catalog, which utilizes a vehicle 
specific debris database to generate propellant debris characteristics. This debris database is a heritage-based/First 
Stage (FS) debris catalog. From this analysis we observed that aborts in the Mach 1 and maximum dynamic pressure 
(Max Q) regions could result in the CM being surrounded by debris for a significant portion of the abort trajectory. 
When the parachutes are deployed, there is a possibility of parachute failure from burn through. Therefore, a 
methodology was developed to calculate the CM parachute temperature profile for a given abort trajectory and a set 
of debris trajectories. 


II. Assumptions/Limitations 

To develop and test the validity of the methodology, several assumptions were made to simplify the problem. 
These assumptions, listed below, represent conservative conditions for the problem. 

• The debris pieces are spherical and radiate as blackbodies. 

• The ambient temperature (T amb ) is 50°F (283 K). 

• The debris burns ( T f ire baii) at a constant 4000°F (2478 K) through the entire abort trajectory. 

• There is no heat transfer between debris, therefore all debris pieces are at the same temperature. 

• Only radiative heating is analyzed; convective cooling is ignored. 

The next set of assumptions is based on Ares 1 and Orion data, which were obtained from the Cx Program. 

• The nylon parachutes have an emissivity (e) of 0.899. 

• The parachutes failure temperature is 150°F (339 K). 

• Debris density is 1811.1 kg/m 3 (used to determine debris size). 

• The sample area of the parachute is, A[ = 1.0 [m 2 ]. 

• The Stefan-Boltzmann constant (o) is 5.669e-9 [W/(m 2 K 4 )]. 

• The Specific Heat of nylon (C) is 1.7 [J/(gK)j. 

• The density of nylon per square area (p N ) is 33.9049 [g/m 2 ] 

III. Methodology 

The foundation of this methodology will be based on the Stefan-Boltzmann Law: 

q = SCXT 4 (Ref. 3) (1) 

The problem is constructed as a three body system, as depicted in Fig. 1. The nylon (body 1) is a section of the 
parachute. The fireball (body 2) is the collection of burning debris. The third body is the portion of the sky that is 
not obscured by the fireball. Our approach was to analyze the heat flux absorbed and expelled by the nylon, by using 
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the Stefan-Boltzmann Law to write two equations: one that represents the heat flux imparted to the nylon from the 
fireball, and one that represents the heat flux radiated out from the nylon to the sky: 


in ~ £ ® fireball 
<lin = £ °{ T nylon 



( 2 ) 

(3) 


Equation (2) assumes all of the radiated heat from the fireball is directed at the nylon. This is not the case, since 
the fireball is a 3D object that radiates in all directions. Equation (3) uses the same assumption, but corresponds to 
the relationship between the nylon and sky. To account for this assumption, Eq. (2) and Eq. (3) will need to be 
multiplied by a view factor. The view factor is. defined as “the fraction of radiation leaving one surface that is 
intercepted by a second surface.” 3 When relating several bodies, the summation rule applies: F 12 +F 13 =l. Applying 
the view factor and summation rule on Eq. (2) and Eq. (3) results in: 

<hn = fireball -KlonYn (4) 

<lout= £(T ( T n y lon - T amb) F n 

k lout = ™{ T nylon ~ T amb )0 ~ F U ) ( 5 ) 


Example: ~ 5° of 
“view” 

(actually is -cone, in 
3D) 


Body 1 - Nylon 
Body 2- Fireball 
Body 3 - Atmosphere 


Sky is remainder; 
»100° of “view” 



Nylon Chute Sample 

Figure 1. Diagram of Thermal Problem 


The system has two equations, Eq. (4) and Eq. (5), and four unknowns: q in , q out , T nylon , and F 12 . Since the goal is 
to create a temperature profile of the nylon, T ny i on will have an initial condition applied based on the situation. This 
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leaves three unknowns, meaning F 12 will have to be solved for outside the heat flux equations. This is accomplished 
by starting with the generic equation for a two body view factor: 


F, 


12 


1 r ("COS#, cos# 2 

~aJ' ~nD 2 


dA ] dA 2 (Ref. 4) 


( 6 ) 


Here Aj is the surface area of body 1, A 2 is the surface area of body 2, 0i is the angle between normal vector of 
A[ and D , 0 2 is the angle between normal vector of A 2 and D , and D is the magnitude of the distance between A] 
and A 2 , An illustration supporting Eq. (6) can be seen in Fig. 2. 



Equation (6) can be integrated to: 


1 


f 




cos 0 1 cos 0 2 

kD 2 


^ 1^2 


( 7 ) 


Now assume that A] is a parachute segment with the normal vector pointed toward A 2 . Only a segment of the 
parachute is analyzed, because even if only one section is degraded, the entire parachute system is considered failed. 
Assume A 2 is the cross section or projected area of a spherical propellant debris piece of radius r. Figure 3 is an 
illustration of F ]2 with these assumptions applied. 



Figure 3. Illustration of F 12 with Assumptions 

Equation (7) will then reduce to: 
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( 8 ) 


Fu = 


cos 0 cos 0 

ZD 2 


-7tr = ■ 


D 1 


The reduction of Eq. (8) implies that the view factor is only dependent on the size of the debris and the distance 
of the debris from the parachutes. F 12 is the view factor of only one debris piece to the parachutes; the view factor 
of all the debris pieces to the parachutes is still needed. I assumed that all the propellant burns at the same constant 
temperature - 4000°F - and that there is no heat transfer between the debris pieces, meaning we can use the 
summation rule - (Fi d i+F ld2 +Fi d3 +Fid 4 + ...) + F [3 = 1, where F [d i is the view factor of body 1 to debris piece 1, Fi^ 
is the view factor of body 1 to debris piece 2, etc - to obtain the F ]2 variable in Eq. (4) and Eq. (5). 

Now that the view factors can be calculated, Eq. (4) and Eq. (5) can be used to calculate the net energy acting on 
the nylon. Multiply Eq. (4) and Eq. (5) by the area of nylon ( A 3 ) to obtain the energy equations for the nylon. The 
net energy on the nylon is then calculated by taking the energy emitted by the nylon and subtracting it from the 
energy absorbed by the nylon: 


i4 rp4 

L fireball ± nylon 




Q»=A T h 

Q„„=ev(T! : „„,-T:,y-F tl )A I 

Qnet Qin Qout 

The net energy can then be used to calculate the change in temperature: 

Qn 


AT = 


c PnA 


(Ref. 4) 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


For each time step, the delta temperature is applied to the previous parachute temperature. The result is a 
parachute temperature profile along the trajectory. An example can be seen in Fig. 5 or Fig. 6. 


IV. Simulations 

The methodology was written into a MATLAB script. The required user inputs are the distance of the debris 
from the CM and the mass of the debris pieces, both at different time steps. The initial temperature of the parachute 
was conservatively set to 70°F (294 K). The debris piece mass is used to calculate the spherical dimensions of each 
piece. The parachute temperature calculator employs the methodology detailed in Section III. The output is a 
parachute temperature profile. Figure 4 displays this simulation process. 


Input; 


Output: 



Figure 4. Simulation Flow Chart 
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V. Results 

Two abort cases were run with the Chute Temperature Calculator, one when the launch vehicle reaches Mach 1 
(—44 sec MET) and one when it reaches Max Q (~60 sec MET). Figure 5 below displays the temperature and 
minimum distance profiles for an abort at Mach 1 . 


Temperature of the Parachutes vs. Time 




Figure 5. Abort at Mach 1 Temperature Profile (a) and Distance of Closest Pieces (b). 


The first temperature increase in Fig. 5 a), from 72 s to 80 s, was caused by a large number of debris pieces 
passing within approximately 200 m to 2000 m of the CM. The maximum temperature seen for the trajectory was 
90°F and caused by the first temperature increase. The circled temperature spikes were caused by one or more 
debris pieces passing within 100 m of the CM. The relation between pieces passing close to the CM and a 
temperature spike can be seen by comparing Fig. 5 (a) and (b) at a selected time - a large downward spike seen in 
the closest debris piece profile correlates to an upwards spike in the temperature profile. The closest debris piece 
from this trajectory was 88 m. For this case, the nylon failure temperature limit of 150°F was not reached. Figure 6 
shows the results for the abort at Max Q. 
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Figure 6. Abort at Max Q Temperature Profile (a) and Distance of Closest Pieces (b). 


The first temperature increase in Fig. 6 (a), from 92 s to 97 s, was caused by a large number of debris pieces 
passing within the approximate range of 200 m to 2000 m. The first temperature spike, at 125 sec, was caused by a 
single debris piece flying within 14.8 m of the CM and only took 1 s to occur. The other smaller spikes were caused 
by a few pieces getting within 100 m of the parachute. For this case, the nylon failure temperature limit of 150°F 
was not reached. 

Even though the Mach 1 and Max Q cases evaluated did not result in a thermal failure of the parachutes, there 
still exists a probability of the parachute temperature reaching the nylon failure limit for other debris field 
configurations. These two cases also demonstrate the risk that large temperature increases can occur in a short 
amount of time. 


VI. Conclusions 

These two cases indicate there is a possibility of the temperature profile reaching the nylon failure temperature, if 
one or two pieces get too close to the CM. Only two cases were run and the results do not represent the full 
spectrum of possible debris field configurations. From these two cases I observed that the possibility exists for large 
temperature increases in a short amount of time, meaning that for a different debris field configuration and dispersed 
abort trajectory, a failure could occur. A set of dispersed abort trajectories with dispersed debris field configurations 
should be analyzed using the developed thermal methodology; however a Monte Carlo of this scale is very 
computationally intensive, and was not performed here due to the limited computational resources available to the 
author. This additional work would be used to calculate the probability of a thermal limit violation for the 
parachutes. 


Acknowledgments 

The author wishes to thank Sean Stapf and Ken Gee for their help with developing the methodology, Mary Jane 
O’Rourke for providing the constants and values used in the simulations and Stephen Miller for reviewing the 
finalized methodology and providing verification of assumptions. 


7 

American Institute of Aeronautics and Astronautics 


References 

1 Chapman, Alan J., Heat Transfer, 4 th Ed., Macmillan, New York, 1984, Chap 11. 

“ Reynolds, William C., and Perkins, Henry C., Engineering Thermodynamics, McGraw-Hill, New York, 1977, pp. 576-580. 

3 Incropera, Frank P., and DeWitt, David P„ Fundamentals of heat and Mass Transfer, 5 th ed., John Wiley & Sons, New 
York, 2002, Chaps 12, 13. 

4 Avallone, Eugene A., and Baumeister, Theodore, III, Marks ’ Standard Handbook for Mechanical Engineers, 10 th Ed., 
McGraw-Hill, New York, 1996, pp. 4-63, 4-64, 4-65. 


American Institute of Aeronautics and Astronautics 



First Stage Solid Propellant 
Multiply Debris Thermal Analys s 


Ben Toleman 


7 / 29/2011 

Debris_Thermal_Analysis.pptx 



USA 

nt: 

1 United Space Alliance 


Outline 


• Introduction 

• Method 

• Simulation 

• Conclusion 


7 / 29/2011 

Debris_Thermal_Analysis.pptx 



Ascent / Descent 
-light Design 


USA 

United Space Alliance 


Introduction 


* A solid rocket stage on a launch vehicle can create a thermal radiation hazard for an 
aborting crew module. 

- For vehicles launched from the Eastern Range at Cape Canaveral, the Air Force requires 
capability from the ground to terminate thrust for launch vehicle on an errant trajectory in 
order to prevent that vehicle from traversing over land and putting the public at risk. 

- Because solid propellant vehicles are incapable of being shut down, they require a Flight 
Termination System (FTS) that, when activated, destroys the vehicle and freezes the 
vehicle’s impact point. 

- Detonation of Constellation’s Ares I launch vehicle results in thousands of debris 
fragments, the majority of which are first stage solid rocket propellant fragments, which 
pose a threat to the aborted crew module (CM). 

- The risk of debris strikes in first stage abort scenarios have been studied in depth by the 
Cx Program. However, a secondary debris concern exists - the thermal effects of burning 
solid rocket propellant on the CM, and specifically on the parachute system. 

* I performed an independent analysis to determine the thermal effects of high 
temperature debris on the CM parachute system. 

* From this analysis I developed a process to assess those thermal. 

- Analyzed the cumulative thermal effects on the parachutes to assess risk of parachute 
failure. 

- Integrated the temperature over time from drogue chute deploy. 
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Assumptions 


• First Stage propellant debris 

- spherical shape 

- radiate as black bodies 

• The nylon chutes have an emissivity of 0.899. 

• The chutes failure temperature is 339 °K (150 °F). 

• Conservation of Energy 

• Debris density is 0.065 lbm/in A 3. 

• Ambient Temperature is 283 °K (50 °F). 

• The debris burns at a constant 4000°F (2478 K) through the entire abort 
trajectory. 

• There is no heat transfer between debris, therefore all debris pieces are at the 
same temperature 
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Method - V ew Factor Illustration 


Example: - 5 deg of 


“view” 


(actually is a -cone 
in 3-D) 


Sky is remainder; 
»100° of “view’ 



Body 1 - Nylon 
Body 2 - Fireball 
Body 3 - Atmosphere 


Nylon Chute Sample 
(~1” square “postage stamp”) 
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Method 


* The foundation of this methodology is based on the Stefan-Boltzmann Law: 

q = soT 4 




The problem is constructed as a three body system, as depicted in the next 
slide. Two equations were written to represent the heat flux imparted to the 
nylon from the fireball, and the heat flux radiated out from the nylon to the 
sky: 



— £<y(T f ireball 



) 



• The two equations above are multiplied by a View Factor to account for the 
percent of radiation transfer. The summation rule, F 12 + F 13 = 1, is applied to 
consolidate variables. 


q in 

Qout 

Qout 


fireball 


T 4 V 7 

± nylon 12 
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Method - Continued 


* F 12 is defined as the “fraction of radiation leaving surface 1 that is intercepted 
by surface 2”. 


= J_ f rcos(6>)cos(6> 2 ) 

A J 


A\ A 2 


tcD 

r 


• Integrated: F u = 


1 


A 


dA { dA 2 

cos(6 , 1 ) cos(# 2 )^ 


V 


71 D 


A x A 2 


J 



For purposes of this analysis, assume: 

- A 1 is the parachute, with the normal vector pointed toward A 2 . 

- A 2 is the cross section or projected area of a spherical propellant debris piece of 
radius r. 


A 2 = 7TT 
01 = 02 = 0 “ 


_ cos(0)cos(0) 2 

12 nD 2 


r‘ 


D 
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Method - Continued 


• The emissive power equations are multiplied by A 1 to get the energy 
equations of the nylon: 

Qin = & (7 1 fireball ~ ^ nylon )^\lA 
Qout = ^ nylon ~ ^ amb ) 0 - _ ^12 )A 

• Then subtract Q in from Q out to get the net energy, Q net . 

Qnet Qin Qout 

• The net energy can then be used to calculate the change in temperature 

4 


A T = 


Qnet 

c PnA 
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Simulation 


* The methodology was written into a MATLAB script. 

* Analyzed the cumulative effects of the changing temperature on the 
parachutes. Integrated the temperature over time from drogue chute deploy. 

* The required user inputs were the distance of the debris from the CM and the 
mass of the debris pieces, both at different time steps. 

* The debris piece mass is used to calculate the spherical dimensions of each 
piece. 

* The methodology detailed in the previous slides was used to create a 
parachute temperature profile. 
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Abort at Mach 1 


Temperature of the Parachutes vs. Time 




Time, s 
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Abort at Max Q 


Temperature of the Parachutes vs. Time 




Time, s 


7 / 29/2011 

Debris_Thermal_Analysis.pptx 


Page 10 



USA 

United Space Alliance 


Conclusions 


* From the two cases run, it was observed that the possibility exists for large 
temperature increases in a short amount of time, meaning that for a different 
debris field configuration and dispersed abort trajectory, a failure could occur. 

* The temperature profile is driven more by a small number of debris pieces 
passing close by (within 100m) than by the general population of debris 
further away. 

* These two cases indicate there is a possibility of the temperature profile 
reaching the nylon failure temperature, if one or two pieces get too close to 
the CM. 

* A set of dispersed abort trajectories with dispersed debris field configurations 
should be analyzed using the developed thermal methodology; however a 
Monte Carlo of this scale is very computationally intensive, and was not 
performed here due to the limited computational resources available to the 
author. This additional work would be used to calculate the probability of a 
thermal limit violation for the parachutes. 
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